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Abstract 


An  alternative  to  Poisson  sampling  called  Sunter  sampling  is  in- 
troduced into  the  forestry  literature.  The  probability  of  selecting  a  sam- 
ple unit  depends  on  its  size  and  the  number  of  previous  units  selected 
as  well  as  on  the  sizes  and  number  of  units  remaining  in  the  popula- 
tion. This  results  in  a  less  variable  sample  size  than  in  Poisson  sam- 
pling. Sunter  sampling  is  more  efficient  than  Poisson  sampling  where 
a  sampling  list  is  available  prior  to  sampling  if  a  slightly  biased 
adjusted  estimator  similar  to  one  in  Poisson  sampling  is  used.  Approx- 
imate true  variances  are  given  for  both  the  unadjusted  and  adjusted 
estimators.  Two  sample-based  variance  approximations  provide  relia- 
ble estimates  of  both  the  true  and  simulation  variance  of  the  adjusted 
estimator.  Sunter  sampling  is  not  yet  a  practical  alternative  when  no 
sampling  list  is  available  but  perhaps  could  be  an  alternative  to,  for 
example,  point-Poisson  sampling. 
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Management  Implications 

A  new  sampling  scheme  is  introduced  that  could  be 
quite  efficient  in  timber  sales  and  other  inventories.  Some 
of  these  inventories  currently  suffer  from  random  sam- 
ple sizes  which  may  result  in  inventories  being  more  cost- 
ly or  less  precise  than  planned  for.  The  study  suggests 
applications  where  this  new  method  may  be  useful. 

Introduction 

Poisson  sampling  was  introduced  by  Hajek  (1957, 
1964)  and  Grosenbaugh  (1964)  as  a  pps  without  replace- 
ment sampling  scheme  when  no  sampling  list  is  avail- 
able. Units  are  selected  exactly  proportional  to  size. 
Because  the  sample  size  is  random,  however,  the  sam- 
ple size  is  frequently  different  from  what  was  intended. 

The  purpose  of  this  note  is  to  discuss  Sunter  sampling, 
an  alternative  for  which  the  sample  size  is  less  variable, 
and  to  compare  it  to  Poisson  sampling  in  applicability 
and  efficiency.  The  objective  is  to  encourage  practition- 
ers to  look  for  practical  applications  for  Sunter  sampling. 
One  possibility  might  be  point-Sunter  sampling  as  an 
alternative  to  point-Poisson  sampling. 

Review  of  Literature 

Sunter  (1986)  gives  an  alternative  to  Poisson  sampling. 
He  starts  with  the  premises  that: 

1.  units  are  selected  with  probabilities,  IIj,  propor- 
tional to  Xj 

2.  the  joint  inclusion  probabilites      of  units  indexed 
i  and  j  are  positive  for  all  i,  j 

3.  nTIj  -       >  0  for  all  i,  j  (i*  j) 

4.  IIj:  is  known  or  can  be  calculated  for  all  i,j. 

Note  that  Poisson  sampling  satisfies  all  these  condi- 
tions except  that  Iljj  =  niITj  for  all  i  ^  j,  so  that  nTIj 
-  Iljj  =  0,  for  all  i*  j. 

As  Sunter  (1986)  indicates,  under  any  set  of  positive 
selection  probabilities,  a  Horvitz-Thompson  estimator 

N 

of  the  population  total  Y  =  L  Yi  is  given  by 

Y  =  .^Yi/IIi  [1] 
where  na  is  the  achieved  sample  size,  with  variance 

VfY)  =.EiYi2(i-ni)/ni-E  ¥^(^11,-11^/114111 

and  unbiased  variance  estimator  (if  Ely  >  0  for  all  i,j) 


v(Y)  ^E^il-U^-'L  Y^U^-U^/U^U^  [2] 
For  fixed  n ,: 


N  N  N 


VfY)  =  £.  Pinj-ni^Yi/ni-Yjiij)2  [3] 

and 

v(t)  =  .^(njnj-ni^Yi/ni-Y/n^/nij  [4] 

Note  that  the  first  premise  arises  from  the  desire  to  make 
the  square  term  in  equation  [3]  small,  which  will  be  true 
if  there  is  proportionality  between  the  variate  value  and 
size  measure.  Satisfying  the  second  and  fourth  premises 
ensures  that  either  equation  [2]  or  equation  [4]  can  be 
computed  as  appropriate.  The  third  premise  implies  that 
the  variance  estimate  in  equation  [4]  cannot  be  negative 
and  indicates  some  stability  in  the  estimate. 

Stuart  (1963)  indicates  that  it  is  advantageous  to  have 
Ili  nj-nj^O  for  the  largest  (Yj/nj-Yj/np2  so  that  there 
are  pps  schemes  that  can  be  more  efficient  than  either 
Poisson  or  Sunter  sampling.  But  it  is  advantageous  to 
be  able  to  compute  all  joint  probabilities  of  selection  and 
to  always  have  positive  variance  estimates.  Sunter  (1977) 
first  proposed  a  list  sequential  sampling  scheme  that 
results  in  fixed  sample  size.  The  disadvantage  of  this 
scheme  is  that  the  second  part  of  the  sample  is  chosen 
by  simple  random  sampling  and  that  the  joint  probabil- 
ity of  selection  is  not  available  between  the  sample  units 
in  the  second  part  and  those  in  the  first  part  of  the 
sample.  Sunter  then  proposes  a  modified  sequential 
sampling  scheme  as  follows: 

a.  Select  a  number  x*  such  that  for  11;  =  ntXj/Z 

N 

ntXi<Xi=  E^Xj  +  x*  (i  =  1,2,. ...N)  [5] 
For  x*  =  0,  ntxn(  < X*(  if  nt  =  1 

where  nt  is  the  target  sample  size.  Sampling  is  propor- 
tional to  x. 

b.  Select  the  i-th  unit  with  conditional  probability 

P(Ui|ni)  =  nixi/Xi*  [6] 

where  Uj  is  the  event  that  unit  i  is  selected,  and  nj 
denotes  the  number  of  units  remaining  to  be  selected 
when  the  i-th  unit  is  reached. 

Sampling  stops  when  nt  units  are  selected  or  when 
the  population  of  units  has  been  exhausted.  Thus,  the 
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achieved  sample  size  na  is  no  greater  than  the  target 
sample  size  nt.  The  above  is  referred  to  as  Sunter  sam- 
pling in  the  rest  of  this  paper. 
Assume,  without  loss  of  generality,  that  the  size  meas- 

N 

ures  x;  are  scaled  such  that   E  x.  =  l,  then  units  are 

i  =  l  1 

selected  with  probability 

n^ntx/X*  (i  =  l,2,...,N) 
and  joint  probabilities 

11^  =  nt(nt-l)Xixj7i 7X*  (j  > i  =  1 ,2 , . . .  ,N-1) 

N 

where  X*  =  X,  *  =  1  +  x*,  X.*  =    E  x;  +  x*, 

j  =  i  ' 

Yl*  =  1/X2*,7i*  =  (l/Xj*+  Jfl-Xi/Xa*)  ... 

(l-xi_1/Xj)(i  =  2,3,...,N-l), 

and  x*  is  a  constant  such  that  ntXj  <  Xj*  for  i  =  1 ,  2 ,  . . . ,  N. 
Note  that  the  expected  sample  size,  ne,  is 

E(na)  =  nt/X*  =ne, 

V(na)  =  .E  1^(1-11^-2  E.  (^0:-^=), 
1  =  1  i<i        '  > 


and 


P(na<nt)<ntx*/X* 


This  sampling  scheme  reduces  the  variance  in  the  sam- 
ple size  na  relative  to  Poisson  sampling  by  the  difference 
of  the  sums  of  product  and  joint  probabilities. 

The  selection  algorithm  will  work  for  any  ordering  of 
the  population,  but  will  be  easier  in  certain  situations. 
An  ideal  ordering  would  be  close  to  the  requirements 
of  Corollary  1.1  in  Sunter  (1986)  which  states:  if 

N 

xN-nt  +  l=xN-nt  +  2  =  ---=xN>   ntxi<Xi  =  .£  xj  with  i  = 

j  =  1 

1,  2,  N-nt  for  some  ordering  of  the  population,  then 
x*  =  0  will  work.  Then  a  sample  of  exactly  nt  units  will 
be  obtained.  This  is  rarely  satisfied  in  practice.  The  next 
best  ordering  may  be  to  arrange  the  population  units  by 
decreasing  order  of  the  size  of  x.  If  a  sampling  list  is 
available,  x*  can  be  computed  easily. 


Analytical  Comparisons 

The  Horvitz-Thompson  estimator  for  Sunter  sampling 

is 

n  _  n  „  Y: 

Yus  =  LaYJU:  =  Ea  — X7nt 
us      i=l    1     1     i=l  Xj  1 

and  can  be  rewritten  as 

Yus  =  X/ne  EaVxi 
us  e.=1 


[7] 


[8] 


where  X  is  the  population  total  of  the  auxiliary  variable  x. 

Analogous  to  the  (biased)  adjusted  estimator,  Ya,  sug- 
gested by  Grosenbaugh  (1967),  a  more  efficient  but  biased 
estimator  should  be 


^as       Yus  **e^a 


The  variance  of  Yus  can  be  found  in  most  sampling 
textbooks.  The  true  variance  of  Yas  is  very  complicated. 
An  approximate  variance  of  Yas  similar  to  the  one  sug- 
gested by  Schreuder  et  al.  (1971)  for  Poisson  sampling  is 


VG(Yas) 


N     X,  Y:X 


i+ 


V(na) 


[9] 


The  most  suitable  situation  for  pps  sampling  without 
replacement  with  the  Horvitz-Thompson  estimator  in 
equation  [1]  is  the  one  suggested  by  Van  Deusen  (1987): 


Yj  =  jfixj  +  ej 
with  Eei  =  0  and  Ee^;  =  J  <t2Xj2  for  i  =  j 


[10] 


Under  this  model 
V(Ya)  =  a2X2/ne 


0  for  i  *  j 


1  +  l/np-  E  x2/X2 
e  i=i  1 


[11] 


where  Ya  is  the  adjusted  estimator  of  Poisson  sampling, 
and  ne  is  the  expected  sample  size.  It  is  shown  in  the  ap- 
pendix that  under  this  same  model 

V(YUS)=  -  (a2  +  i32)X*X+  —  /J2X*X  E  7l*  x,  [12] 


i  =  l 


ru-1 


-  -^-{o2  +  I32)X*  X^xf-p  X2 
and  an  approximate  variance  for  Yas  is 


V(Yas)  =  a2X2/nf 


1  N 
 nt(nt-l)  (X2/X*-E  x2/X7Tl*) 


i=i 


[13] 


The  comparison  does  not  lead  to  a  simple  analytical  com- 
parison of  which  method  (Poisson  or  Sunter)  has  the 
smallest  variance  because  the  joint  probabilities  of  selec- 
tion under  Sunter  sampling  are  analytically  very  com- 
plex and  no  simplifying  assumption  is  obvious.  The  two 
methods  are  compared  based  on  their  approximate  true 
variances  and  also  by  a  simulation  study  for  some  popu- 
lations. The  approximate  true  variance  of  Yas  is  given  in 
the  appendix. 

Two  sample-based  variance  approximations  suggest 
themselves  for  V(Yas).  The  first  one  is  the  Yates-Grundy 
variance  estimate  which  assumes  a  fixed  sample  size,  i.e., 
v(YYg)  =  v(Y)  in  equation  [4]  and  the  second  one  is  with 
the  replacement  variance  estimator  corrected  for  the  ran- 
dom sample  size 


X2     na  Y: 

vr(Yas)  =           Ea(~^  -  E 

ne  i  =  i  Xj     i  =  i  neXj 


)2/(ne-l)  [14] 


Populations  and  Simulations 

To  develop  an  understanding  of  gains  and  losses  in 
efficiency  and  bias  due  to  the  validity  of  error  structure 
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for  postulated  linear  models,  the  loblolly  pine  popula- 
tion described  in  Schreuder  and  Thomas  (1985)  was  used 
as  a  starting  point.  This  population  consists  of  1,801  lob- 
lolly pine  trees  with  cubic  foot  volume  as  the  dependent 
variable  (y)  and  tree  basal  diameter  squared  times  total 
height  (x)  as  the  auxiliary  variable.  This  data  set  is  basi- 
cally a  subset  of  a  larger  data  set  of  5,500  trees  described 
in  McClure  et  al.  (1983),  selected  such  that  the  diameter 
distribution  of  these  trees  is  similar  to  that  for  loblolly 
pine  in  the  Piedmont  of  the  Southeast,  a  highly  skewed 
distribution  with  many  small  and  few  large  trees. 

From  this  data  set  (N  =  1,801)  the  following  parameters 
were  computed.  The  mean-of-ratios  regression  parameter 
is 

0mr  =  JiYi/Nxi 

with  variance  cx2mr=  L  [{Yr0mixiyx]\£/N 
The  ratio-of-means  regression  parameter  is 

N  N 
N 

with  variance  a2Tm  =  L  [(Yj-^Xj)  /  Vx[]2/N 
The  simple  linear  regression  parameter  is: 

(3T= J^v?)  (Xi-xyJjXi-xp 

with  intercept  parameter  ar  =  Y-jSr  X.  Variances  used 
with  these  two  parameters  were 

N  r 

a2rl  =  E  [(Yj  -  ar  -  ^J/xJ  2/N  and 
ff2r2  =  Jj(Yi-ar-^i)/VxT]  2/N 

Then  populations  1-10  were  generated  in  pairs  (repli- 
cates) using  the  same  x's  as  in  the  original  population 
and  with  the  following  Y's. 

For  populations  1  and  2: 
Yi  -  ej  =  /3mr  Xj  where  et  ~  N  (0,  a2mT  xf) 

For  populations  3  and  4: 
Yi  -  ej  =  ^  Xj  where  e;  ~  N  (0,  o2Tm  Xj) 

For  populations  5  and  6: 
Yj  -  e;  =  ar  +  |8r  xj  where  ej  ~  N  (0,  a2rl  x;2) 

For  populations  7  and  8: 
Yj  -  e;  =  ar  +  j8r  xi  where  ej  ~  N  (0,  a2r2  Xj) 

and  for  populations  9  and  10: 
Yj  -  ej  =  ar  +  (8r  Xj  where  ej  ~  N  (0,  a2r2  Xj21) 

To  save  computer  time  the  same  random  numbers  were 
used  for  populations  1,  3,  5,  7,  and  9  and  similarly  for 
populations  2,  4,  6,  8,  and  10.  Population  11  is  the  origi- 
nal loblolly  pine  population.  All  populations  are  the 
same  size.  Parameter  values  are  (3mT  =  0.001836,  (3Tm  = 
0.002016,  aT  =  -1.146,  and  j8r  =  0.002209. 


Another  population  data  set,  the  Pine  Creek  diameter 
data  set,  for  which  2,748  tree  diameters  (D)  were  meas- 
ured in  1967  and  1979,  was  also  used.  The  size  of  the 
intercept  is  close  to  the  size  of  the  slope  for  the  relation- 
ship between  D2  in  1979  and  1967.  The  plotted  relation- 
ship seems  to  indicate  a  very  strong  linear  relationship 
with  variability  increasing  with  increase  in  diameter 
class  size. 

For  the  Pine  Creek  data  set,  populations  1-10  were 
generated  as  above  for  the  loblolly  pine  but  with 
parameter  values  /3mr  =  1.323,  0^  =  1.272,  ar  =  2.75, 
and  pT  =  1.110. 

As  with  the  loblolly  pine,  the  same  random  numbers 
were  used  for  populations  1,  3,  5,  7,  and  9  and  for  2,  4, 
6,  8,  and  10,  respectively.  Population  11  is  the  original 
Pine  Creek  diameter  population.  All  11  populations  are 
the  same  size. 

Using  Poisson  sampling  and  Sunter  sampling,  2,000 
samples  of  intended  size  20  were  drawn  from  each  popu- 
lation. In  Sunter  sampling,  the  order  of  visiting  the  units 
is  very  important.  First  units  were  arranged  in  decreas- 
ing size  of  the  auxiliary  variable  with  x*  =0.12%  of  X  for 
the  loblolly  pine  and  x*  =0.37%  of  X  for  the  Pine  Creek. 
This  is  a  good  strategy  for  minimizing  the  value  of  x* 
so  as  to  guarantee  TIj  =  nxj/X*  <1  for  all  i  (i=l,  N). 
The  expected  sample  size  is  19.98  for  the  loblolly  pine 
(x*  =0.12%X),  and  19.93  for  the  Pine  Creek  (x*  =0.37%X). 

To  evaluate  other  orderings  of  the  population,  the 
population  units  were  randomly  ordered  twice.  For  these 
orderings  x*  =4.27%X  and  x*  =5.23%X  were  required  for 
these  two  situations  for  the  loblolly  pine  data  set  and 
x*  =  1.05%X  and  x*  =  1.26%X  for  the  Pine  Creek  data 
sets.  The  expected  sample  sizes  are  19.15  (x*  =4.27%X) 
and  19.02  (x  =5.23%X),  respectively,  for  the  two  loblolly 
pine  data  sets,  19.79  (x*=1.05%X)  and  19.75 
(x*  =1.26%X),  respectively,  for  the  two  Pine  Creek  data 
sets. 


Computational  and  Simulation  Results 

In  all  cases  in  the  simulation  study,  the  estimation 
biases  of  Yus,  Yas,  and  Ya  were  negligible  and  are  not 
shown.  In  table  1,  the  expected  variance  under  the  model 
YUS(V(YUS))  is  shown  in  equation  [12],  the  approximate 
true  variance  for  Yas(V(Yas))  in  equation  [13],  the  approx- 
imate true  variance  for  Ya(V(Ya))  in  equation  [11],  simu- 
lation variances  for  Yas  and  Ya  based  on  2,000 
simulations  and  average  sample-based  variance  estimates 
for  Yas  based  on  fixed-sample  size  (v(YYg))  in  equation 
[14],  and  based  on  with-replacement  sampling  (vG(Yas)) 
in  equation  [14]  for  populations  1  and  2  of  each  data  set 
(only  these  two  populations  were  generated  using  equa- 
tion [4]).  Based  on  the  theoretical  variances  it  is  clear 
that  Yas  is  more  efficient  than  Yu_  for  Sunter  sampling 
and  also  more  efficient  than  Ya  for  Poisson  sampling. 
The  simulation  variances  confirm  that  Yas  for  Sunter 
sampling  is  more  efficient  than  Ya  for  Poisson  sampling. 

Either  v(YYg)  or  vG(Yas)  seem  to  give  reliable  variance 
estimates  of  the  true  variance  of  Yas  (expressed  either  by 
the  approximate  variance  V(Yas)  or  the  simulation  vari- 
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Table  1.— Comparison  of  theoretical,  iterated,  and  average  sample  variances  for  Sunter  and  Poisson 
sampling  for  populations  1  and  2  of  both  data  sets  with  x  sorted  in  decreasing  order. 


Data  set  Method  x* 


Loblolly               Sunter  0.12%X 

Sunter  0.12%X 

Poisson   

Sunter  0.12%X 

Poisson   

Sunter  0.12%X 

Sunter  0.12%X 

Pine  Creek          Sunter  0.37%X 

Sunter  0.37%X 

Poisson   

Sunter  0.37%X 

Poisson   

Sunter  0.37%X 

Sunter  0.37%X 


Population 

Estimator  Variance  1  2 


V(Yus)/106 

0.89 

0.85 

Yas 

V(Yas)/106 

0.79 

0.78 

Ya 

V(Ya)/10b 

0.85 

0.80 

Yas 

Iterated 
variance/10b 

0.74 

0.78 

Ya 

Iterated 
vanance/10b 

0.82 

0.79 

Ya 

v(YG)/10b 

0.78 

0.79 

Ya 

VG(Yas)/1° 

0.80 

0.76 

Yas 

V(YUS)/107 

1.06 

1.07 

Yas 

V(Yas)/107 

0.82 

0.81 

Ya 

V(Ya)/107 

0.89 

0.85 

Yas 

Iterated 
variance/107 

0.85 

0.85 

Ya 

Iterated 
variance/107 

0.88 

0.83 

Ya 

v(YG)/107 

0.82 

0.89 

Ya 

vG(Yas)/107 

0.85 

0.82 

ance).  Neither  sample-based  variance  estimator  is  clear- 
ly better. 

Table  2  shows  a  very  consistent  though  small  advan- 
tage in  efficiency  for  Yas  (Sunter  sampling)  relative  to  Ya 
(Poisson  sampling)  for  almost  all  of  the  11  populations 
in  each  data  set.  This  is  true  even  though  the  expected 
sample  size  for  Sunter  sampling  is  somewhat  smaller  es- 
pecially for  the  x*  values  much  larger  than  0.  There  is 
little  difference  in  the  results  for  the  three  different  x* 
values  even  though  the  expected  sample  size  for  Sunter 
sampling  is  less  with  increasing  x* .  There  is  no  obvious 
change  in  gain  in  Sunter  sampling  with  randomly  or- 
dered populations  over  Poisson  sampling  for  different 
variance  functions  or  presence  or  absence  of  an  intercept 


in  the  y-x  relationship.  Estimator  Yus  of  Sunter  sampling 
could  be  very  inefficient  especially  when  the  x*  values 
are  much  larger  than  0.  As  could  be  expected,  there  seems 
to  be  some  correlation  in  the  results  of  populations  1,  3, 
5,  7,  and  9  and  those  of  populations  2,  4,  6,  8,  and  10. 

Practical  Considerations 

If  a  sampling  list  is  available  prior  to  sampling,  there 
is  no  question  that  Yas  for  Sunter  sampling  should  be 
preferred  over  Ya  for  Poisson  sampling.  By  arranging  the 
population  in  decreasing  order  of  x,  x*  can  usually  be 
taken  to  be  close  to  zero  so  that  the  attained  sample  size 


Table  2.— Ratio  of  square  root,  iterated  mean  square  errors  of  Poisson  over  Sunter  sampling  ratio  estimators  for  the  loblolly  pine-based  and  Pine  Creek 
diameter-based  populations  for  x  sorted  from  largest  to  smallest,  and  for  two  randomly  arranged  sequences  with  2,000  iterations. 


Population 

a  =  0  a  =  0 

Data  Set  x*  Estimator        1  23456789  10  11 


Loblolly 


Pine  Creek 


0.12%X 
4.27%X 
5.23o/oX 
0.37%X 
1.05%X 
1.26%X 


Yas 
us 
as 
us 
as 
us 
as 

r: 
h 


1.00 

1.03 

1.02 

1.03 

1.02 

1.02 

1.03 

1.02 

1.02 

1.02 

1.01 

0.98 

1.02 

0.99 

1.01 

1.01 

1.01 

1.01 

1.01 

1.01 

1.01 

0.99 

0.99 

1.03 

0.97 

1.03 

1.01 

1.01 

0.99 

1.02 

1.01 

1.01 

1.01 

0.69 

0.73 

0.60 

0.63 

0.72 

0.74 

0.67 

0.69 

0.74 

0.76 

0.63 

1.01 

1.01 

1.02 

1.01 

1.04 

1.01 

1.04 

1.00 

1.03 

1.01 

0.98 

0.70 

0.67 

0.61 

0.58 

0.72 

0.70 

0.67 

0.63 

0.74 

0.70 

0.58 

1.00 

1.03 

1.00 

1.02 

1.01 

1.03 

1.00 

1.02 

1.01 

1.03 

1.03 

0.89 

0.90 

0.88 

0.87 

0.90 

0.91 

0.88 

0.88 

0.90 

0.91 

0.87 

1.03 

0.99 

1.03 

0.99 

1.03 

0.99 

1.03 

0.99 

1.03 

0.99 

1.03 

0.75 

0.70 

0.72 

0.67 

0.76 

0.71 

0.73 

0.70 

0.77 

0.72 

0.72 

1.07 

1.03 

1.06 

1.02 

1.09 

1.04 

1.08 

1.02 

1.09 

1.03 

1.04 

0.73 

0.72 

0.70 

0.69 

0.76 

0.74 

0.72 

0.70 

0.77 

0.74 

0.67 
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will  be  identical  to  or  nearly  identical  to  the  desired  sam- 
ple size.  The  simulation  results  show  that  Yas  for  Sunter 
sampling  will  be  more  efficient. 

When  no  sampling  list  is  available,  for  example,  in  tim- 
ber sales  for  which  Poisson  sampling  was  originally  sug- 
gested in  forestry  (Grosenbaugh  1964),  it  is  unlikely  that 
Sunter  sampling  is  a  practical  alternative.  Imple- 
mentation of  this  method  requires  the  ability  to  compute 
the  probability  of  selecting  a  unit  which  depends  on  the 
number  selected  and  the  number  and  sizes  of  units  re- 
maining to  be  selected  in  the  population,  as  indicated 
by  equations  [5]  and  [6].  This  would  likely  be  too  cum- 
bersome in  practice,  particularly  in  light  of  the  small  gain 
in  efficiency  this  seems  to  yield  as  indicated  in  table  1 
for  the  x*  values  considered. 

Perhaps  point-Sunter  sampling  might  be  a  practical  al- 
ternative to  point-Poisson  sampling.  Point-Poisson  sam- 
pling (Grosenbaugh  1971,  1979;  Schreuder  et  al.  1984)  is 
a  highly  efficient  sampling  strategy  where  trees  at  each 
point  are  first  selected  proportional  to  basal  area  (point 
sampling)  and  then  a  subsample  of  trees  is  selected 
proportional  to  tree  height  by  Poisson  sampling  for 
volume  measurement.  In  some  cases,  it  may  be  desira- 
ble to  select  one  or  more  sample  trees  per  sample  point. 
This  is  not  guaranteed  by  point-Poisson  sampling  but 
might  be  guaranteed  by  implementation  of  some  point- 
Sunter  sampling  scheme  (where  the  subsample  of  trees 
at  each  point  is  selected  by  Sunter  sampling).  For  each 
sample  point,  a  list  of  trees  could  be  obtained  so  that  Sun- 
ter sampling  could  be  implemented  at  those  points  where 
sample  trees  are  desired.  This  would  probably  require 
a  programmed  computer  to  indicate  sample  trees.  The 
results  of  the  simulations  offer  some  encouragement  that 
this  may  also  be  more  efficient  than  point-Poisson 
sampling. 
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Appendix 


It  is  known  that  the  variance  for  the  Horvitz -Thompson 
estimator  is 


N  1-FL         N  N  IL;-rL  n 
V(YUS)  =L  -j^Yf+L  L 


Y-  Y- 


j 

N    Y2        N  N  I!.-; 
=  ?  "FT  +^?TTTTYiYi-(?Yi)2 


i  rii     igfcj  n4  n 


Since 


n^nj  Xj/X*  and  ni:=nt(nt-l)  xi  Xj  7-7X* 

for  Sunter's  sampling,  we  have 

N   Y2X*        N  N  nt-l  N 

V(YUS)=  E  —         +  E  E  Tl*  X^YpfE  Y^ 

i     ntXj        i*j     nt  '  i 

N   Y2X*      nt-l      N  N 

=  ?  -7—  +  -J-X*(E  TiYjY-E  TfYf)-Y2 
1     ntXj         nt        1  1 

where  Y=  E  Y:. 

i 

We  assume  that  Xj  and  7i  are  constants,  and  under 
equation  [10] 


E(Yi)  =  ^xi       E(Y2)  =  (<x2  +  /32)x2 


Hence 


(<72  +  i32)  x?  X*       nt-l  , 

v(Yus)=  ?-  ^r2 —  +  -Vx  /3X  e  7i/3  Xi 


 ^ — X*  E  7*((j2  +  /32)  x2-/32  X2 

(a2  +  02)  X*X     nt-l    ,  , 

=  e{  +  ^2xx?^  xi 


iu-i 


-  — —  (a2  +  02)  X*  E  7i*  x?-/32  X2 


To  derive  the  variance  of  Yas,  the  conditional  variance 
is  used  as  follows: 

V(Yas)=V[E(Yas|na)]+E[V(Yas|na)] 


K    -  1 

V 

^E(Yus|na) 
■  a 

+  E 

-2  V(Yus|nJ 

■  a 

Hence  the  relationship  Yas  =  Yus  ne/na  is  used.  When  na 
is  fixed,  we  have 


E(Yus|na)=Y 


and 


V(Yus|nJ=E  E  (rqiij-iiij)  (Y/ni-Y/n^ 


The  above  two  quantities  on  the  right  hand  side  of  the 
equations  do  not  depend  on  na.  So  by  Taylor  series  ap- 
proximation,we  have 


and 


where 


Hence 


V(-^  =n2v(-)  =n2-^- 
V      e   V      e  n| 


E(4)=n2(-2   +3  -y^), 
n2  n2  n| 


v(na)=E  iii  (i-iy  -2  e  e  (n^-n^ 


V(na)  v(na) 

v(Yas)^Y2  — r-  +[i+3  -f )  e  E^nj-nip 


ne       1  <) 


(Yi/ni-YjH)2 
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Sunter  sampling  is  an  unequal  probability  sampling  procedure 
where  the  probability  of  selecting  a  sample  unit  depends  on  its  size, 
the  number  of  previous  units  selected,  and  the  sizes  and  number  of 
units  remaining  in  the  population.  This  method  is  compared  to  Pois- 
son sampling. 
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